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We derive an operator identity which relates tight-binding Hamiltonians with arbitrary hopping 
on the Bethe lattice to the Hamiltonian with nearest-neighbor liopping. This provides an exact 
expression for the density of states (DOS) of a non-interacting quantum-mechanical particle for 
any hopping. We present analytic results for the DOS corresponding to hopping between nearest 
and next-nearest neighbors, and also for exponentially decreasing hopping amplitudes. Conversely 
it is possible to construct a hopping Hamiltonian on the Bethe lattice for any given DOS. These 
methods are based only on the so-called distance regularity of the infinite Bethe lattice, and not on 
the absence of loops. Results are also obtained for the triangular Husimi cactus, a recursive lattice 
with loops. Furthermore we derive the exact self-consistency equations arising in the context of 
dynamical mean-field theory, which serve as a starting point for studies of Hubbard-type models 
with frustration. 

PACS numbers: 71.10.Fd, 71.27.-fa, 05.50.-fq, 02.10.Ox 
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I. INTRODUCTION 

The Bethe lattice is an infinite graph, where any two 
points are connected by a single path and each vertex 
has the same number of branches Z, as shown in Fig.^ 
for Z = 4. The name "Bethe lattice" originates from the 
fact that Bcthc's approximation for the Ising model is ex- 
act on this latticei^ A finite portion of the Bethe lattice 
is called Caylcy tree. The latter has a peculiar thermo- 
dynamic limit due to its large surface fSii^ whereas the 
infinite Bethe lattice has no surface, all its lattice sites 
being located inside the infinite tree. 

Strictly speaking the Bethe lattice is a pseudolattice 
because it does not possess the usual point and trans- 
lational symmetries of crystal (Bravais) lattices. Nev- 
ertheless, it plays an important role in statistical and 
condensed-matter physics because some problems involv- 
ing disorder and/or interactions can be solved exactly 
when defined on a Bethe lattice, e.g., Ising modelsfi*2i& 
percolation^iSiS or Anderson localizationiiSiiiiiSii^ Such 
exact solutions on the Bethe lattice for Z < oo 
sometimeSji^ but not always fi2*ii*i2ii^ have mean-field 
character. Furthermore, it was argued that mean- 
field theories are more reliable if derived on a Bethe 
latticciii On the other hand, the Bethe lattice may ac- 
tually serve as a model for the electronic structure of 
amorphous solids, as proposed by Weaire and Thorpe^^ 
e.g., for hydrogenated amorphous siliconii^ Laughlin and 
Joannopoulosii used the Bethe lattice to describe lat- 
tice vibrations in amorphous silica, and recently this ap- 
proach was also applied to phonon transport through 
silica-coated nanowircsi^ 

There are two special properties that make the Bethe 
lattice particularly suited for theoretical investigations. 



One is its self-similar structure which may lead to re- 
cursive solutions. The other is the absence of closed 
loops which restricts interference effects of quantum- 
mechanical particles in the case of nearest-neighbor (NN) 
coupling. The situation is different if also longer-range 
hopping processes or interactions are allowed, e.g., be- 
tween next-nearest neighbors (NNN). For example, the 
frustration introduced by NNN hopping typically sup- 
presses antiferromagnetism in the half-filled Hubbard 
model at weak coupling. 

In this paper we consider tight-binding Hamiltonians 
describing hopping of a single quantum particle, paying 
special attention the case where the hopping has both 
NN and NNN contributions (with respective amplitudes 
ti, 12), and to the limit Z ^ 00. The derivation of 
the spectrum for NN hopping has a long history involv- 
ing many different methods li^iSSiSiiSSiS^ However, these 
methods are not immediately useful for longer-range hop- 
ping. Here we develop a method, based on the algebraic 
properties of so-called distance-regular graphs, which can 
effectively treat arbitrary-range hopping on the Bethe 
lattice and certain other lattices. In particular, we de- 
rive an operator identity [Eq. relating tight-binding 
Hamiltonians H^, which describe hopping between sites 
that are d NN steps apart, to Hi. the Hamiltonian for 
NN hopping. This general result has several applications. 
For comparison with earlier methods we first note that 
the nonlocal Green function for NN hopping can be ob- 
tained via a rather short route. We proceed to derive the 
exact density of states (DOS) for arbitrary hopping and 
discuss in detail the case of t\-t2 hopping and exponen- 
tially decreasing hopping. The inverse problem, i.e., the 
construction of a tight-binding Hamiltonian on the Bethe 
lattice corresponding to a given DOS, is also solved. 
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An important limit for any lattice is that of infinite 
coordination number, Z —^ oo, since it always leads to 
a mean-field theory of some sorti^ It is well known that 
for fermionic lattice models the hopping matrix elements 
must then be properly rescaled, e.g., ti ~ for 
NN hopping.Si In this limit dynamical mean-field the- 
ory (DMFT>2SiSSi2LSiS2i2a becomes exact, yielding self- 
energies that are local in space. In particular the Hub- 
bard model may be mapped onto a single-impurity prob- 
lem with self-consistency conditioni^S Here the Bethe 
lattice leads to further simplifications, e.g., in the so- 
lution of the self-consistency equations, partly due to 
the resulting se mielliptical d ensity of states for NN hop- 
ping, p{e) = 2y/l- {e/D)^/iTTD), where D is the half- 
bandwidthiSii^ Interestingly the algebraic band edges of 
this model DOS resemble those of three-dimensional sys- 
tems. DMFT with this DOS has been instrumental in 
clarifying the phase diagram of the Hubbard model, in 
particular concerning the Mott transition from a para- 
magnetic metal to a paramagnetic insulator at half- 

f^^-^^^,2R.2fi.27.28.23..Sn..^1..^2..S3.34..T-S 

It should be noted that the Mott transition in the Hub- 
bard model with NN hopping is usually hidden by an 
antifcrromagnetic low-temperature phase ^SLSSi^ whereas 
for transition metal oxides such as V2O3 the Mott tran- 
sition line extends beyond the antiferromagnetic phase 
boundary^ In order to describe this experimental phase 
diagram at least qualitatively in terms of a Hubbard 
model, the strong tendency towards antiferromagnetism 
must be reduced, e.g., by including the ever-present NNN 
hopping ^2- This was demonstrated within a DMFT 
setup in which the DOS remains semiclliptic, leading to 
a suppression of antiferromagnetism without modifica- 
tion of the paramagnetic phase ^Si^SiSi^ the results of 
this setup are valid for random hopping on the Bethe 
latticeiS2i2& In the present work we consider standard 
tight-binding hopping without randomness; we find in 
particular that for ti-t2 hopping on the Bethe lattice the 
DOS is no longer semiclliptic but becomes asymmetric. 
We also evaluate the exact DMFT self-consistency equa- 
tions for arbitrary hopping on the Bethe lattice, includ- 
ing phases with broken sublattice symmetry. Our results 
for ti-t2 hopping [Eq. H65() ] differ from the corresponding 
equations employed in Rcfs. i2S,36..38.3!i therefore the 
latter only apply to random hopping. 

The paper is structured as follows: In Sec. |H]we con- 
sider topological aspects of tight-binding Hamiltonians. 
In Sec. mil a general operator identity for the Bethe lat- 
tice is derived, which is used in Scc. lIVI to obtain the DOS 
for various hopping ranges; the Appcndixcontains similar 
results for the triangular Husimi cactus. Hopping ampli- 
tudes for a given DOS arc constructed in Sec. IIVDI In 
Secl^the DMFT self-consistency equations are derived. 
Our results are discussed in Sec. IVII Throughout the pa- 
per we make contact with particular results for the DOS, 
Green functions, and DMFT self-consistency equations 
which were previously obtained by other methods 
a detailed comparison will be discussed in a separate 



publicationisa 

II. TIGHT-BINDING HAMILTONIANS AND 
TOPOLOGICAL LATTICE PROPERTIES 

In this section we discuss relations among tight-binding 
Hamiltonians for lattices that belong to a certain class of 
graphs. In general a graph G is defined as a set of vertices 
V = {j, j, fc, ...} representing nodes or sites of a network, 
and a set of links between these vertices. In particu- 
lar, any Bravais lattice is a graph with vertices and links 
corresponding to lattice sites and bonds between near- 
est neighbors, respectively. For a given lattice one dis- 
tinguishes between the metric and topological distance 
between two vertices i and j. The metric distance is de- 
termined by the metric properties of the space in which 
a graph is embedded. On the other hand, the topological 
distance between sites i and j, denoted by dij hereafter, 
is the smallest number of links joining i and j. In this pa- 
per we only use the topological properties and distances 
of the Bethe lattice. Note that the Bethe lattice can 
be embedded in a hyperbolic (Lobachevsky) space with 
metric properties different from Euclidian spaced 

We now consider hopping Hamiltonians on an arbi- 
trary, infinite graph; for simplicity we assume at most 
one link between two sites and no loops of length one. 
In terms of the quantum-mechanical single-particle op- 
erator which removes a particle from site j and 
recreates it at site i, the general tight-binding hopping 
Hamiltonian can be written as 

H= J2 Uj\^)U\=J2*'iH,, (la) 

ijeV d>0 

where 

H,= Y.\^){j\ (lb) 

hjev 

dij—d 

describes hopping between sites i and j separated by 
topological distance d, i.e., tij = t^. . By definition all 
nonzero matrix elements of Hi are equal to 1, i.e., 

(H.h = m.\j) = [l ^.t^ii; (2) 

and Ho ~ J^i = 1 is the identity. 

The topological properties of a graph are completely 
described by the nearest-neighbor hopping Hamiltonian 
Hi. This concept is used in graph theory where Hi is 
called adjacency matrixm^ Furthermore, a simple inter- 
pretation can be given to the matrix elements of higher 
powers of the hopping Hamiltonian Hi, i.e., (i?")y is the 
number of paths connecting sites i and j in n NN stepsm^ 
Explicitly, one has 

(Hl)rj E (^l)^fci (^1)'=!'=^ • • • ' (3) 

fei,...fc„_iev 
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sponding matrix elements 



FIG. 1: Part of the Bethe lattice with coordination number 
Z = A. Any two sites are connected by a unique shortest path 
of bonds. Starting from the site marked by the open circle, 
horizontally shaded circles can be reached by one lattice step 
(NN), vertically shaded circles by two lattice steps (NNN), 
and doubly shaded circles by three lattice steps. Note that 
the lattice is infinite and that all sites are equivalent; the 
shading appears only for visualization of hopping processes. 



where a term in the sum on the right hand side is equal to 
unity if the string of indices {ikik2 ■ ■ ■ k„-ij) represents 
a path joining i and j, and is zero otherwise. Hence, 
{H")^j is equal to the number of paths connecting i and 
j with n NN steps. 

On the Bcthc lattice the quantity turns out to 

be a function only of n and of the topological distance dij , 
while the specific positions of i and j are unimportant. 
Graphs with this property are called distance regular in 
graph theory4^ The triangular Husimi cactus, which is a 
set of triangles connected by vertices (see the Appendix) , 
is also distance-regular, and hence this property does not 
depend on the absence of loops. However, the situation 
is different for periodic (Bravais) lattices. For example, 
given a site i on the two-dimensional square lattice, there 
are two different positions for a next-nearest-neighbor 
site (dij = 2): one along an axis, the other on the di- 
agonal of a plaquette. In the former case there is only 
one path of length two between i and j, whereas in the 
latter case there are two paths. Thus for crystal lattices 
not only topological distances but also specific positions 
play a role in determining iJ". 

For distance-regular graphs Hd can be written in terms 
of powers of Him^ We show this by first proving the 
inverse relation, i.e.. 



(4) 



d=0 



where alf"^ is the number of paths with n NN steps be- 
tween sites separated by a topological distance d. This 
relation can be easily verified by calculating the corrc- 



(5) 



d=0 



which 



Using {Hd)ij = Sd,d,j one finds (i/Dy 
indeed agrees with Eq. , hence proving . Since 
^ for all 71, and a^n'^ = if d > n, this (triangular) 
system of equations can always be inverted, yielding 



(6) 



Therefore for distance-regular graphs the hopping Hamil- 
tonians Hd are given by polynomials in Hi of order d. For 
example, for a Bethe lattice with coordination number Z 
the first few equations read 



(i7i)2 = Z1 + H2, 
(i?i)3 = {2Z - l)Hi - 
[Hif = Z{2Z - 1)1 



- {3Z - 2)H2 



Ha 



(7a) 
(7b) 
(7c) 



By contrast, such relations do not exist for all path 
lengths n and topological distances d on periodic lattices, 
although low-order relations may be found for certain 
graphs. For example for the honeycomb lattice, which 
is similar to a Bethe lattice with Z = 'i when taking at 
most two NN steps, Hf can be related to Hi and H2 as 
in Eq. (fTafl . 



III. OPERATOR IDENTITIES FOR BETHE 
LATTICES 

We now determine the coefficients in Eqs. (@J and Q 
for the Bethe lattice. These equations can be summarized 
as an operator identity involving the hopping Hamilto- 
nians Hd- As an application the Green function for NN 
hopping is obtained. 



A. Recursive relations and generating function 

We consider a Bethe lattice with coordination number 
Z > 2, i.e., branching ratio K = Z — 1 > 1. In this case 
the coefficients aif^ in are obtained from a simple 
recursion, starting from aQ^' = a\^^ = 1. Namely, each 
path of length n > 1 joining a site i to a different site 
j (with d = dij > 1) is composed of a path that joins i 
to some nearest neighbor of j within n — 1 steps and one 
final step to j. Since for the last step there is only one 
possibility, alf'' is given by the number of possible paths 
with (71 — 1) steps between i and the nearest neighbors 
of j. Of the Z nearest- neighbor sites of j, Z ~ 1 are 
separated from i by a distance d+1 and one by distance 
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d — 1. Together with a similar argument for d = we 
thus have the recurrence relations 







(8a) 






(8b) 




= 1, 


(8c) 



where the last equation is due to the treelike structure of 
the Bethc lattice. Since on the Bcthe lattice it is impos- 
sible to return to the same site within an odd number of 
steps, we note that a^"^ vanishes if n + d is odd. 

These recursion relations for the Bethe lattice can 
be solved in closed form. One first considers the 
generating functions Fd{u) = X^^o which ap- 
pear when summing over Eq. Q), i.e., [1 — uHi]~^ 
= '^'^Q HciFci{u). They obey the recursion relations 
Fd{u) = KuFd+i{u) + uFd-i{u) for d > 1 and Fq{u) 
= 1 + uZFi{u). These equations may be solved by 
the ansata^ Fd{u) = f{u)[ug{u)]'^, which yields g{u) 
= 2/[l + VI -4:Ku^] and/(M) = [1- Zu'^g{u)]-\ where 
the sign in front of the square root has been chosen so 
that Fd{u) = u'^ {1 + 0{u)) is satisfied. By setting x = 
ug{u) and rearranging the terms we finally obtain the 
remarkable operator identity 



1 - x'' 



1 — xHi + Kx^ 



E 

d=0 



Hdx''. 



(9) 



From this formal power series the coefficients in Eqs. (0} 
and l|njl can be extracted in closed form. After some 
algebra we find that the nonzero coefficient a^f' and A^f' 
are given by 



(n-2s) 



E 

r=0 



■2s) 



-irE 



n 
r-l 



n — s — r 
s — r 



(10) 



(11) 



for < 2s < 71. By induction with respect to n it is also 
straightforward to verify that Eq. (|10|) is the solution of 
the recursion relations © . 

Note that in the limit of infinite connectivity, K —* 
oo (or Z oo), one must scale24 Hd by Z^^^^^ where 
Zd is the number of sites that a given site is connected 
to by Hd- To leading order one has Zd ^ Z"^ ^ K"^; 
thus we introduce scaled operators, Hd '■= Hd/K'^^^. By 
computing coefficients we then obtain from Eq. © the 
relation 



Hd = UdiHi/2) 



UniHl/2) = 



1 - ^do - Sdi 
K 

L"/2J ^ 



Ud-2{Hi/2), 



s=0 



(12a) 
(12b) 



which are valid for any K > \ and reduce to Hd = 
Ud{Hi/2) in the limit K oo. Here Un{x) are the 



Chebyshev polynomials of the second kind(^ with gener- 
ating function [1 — 2xt -\- = X]^!°=o Un{x)t'^ ■ It seems 
that the operator identities ((HJ and H12|l might also be 
useful in other contexts involving the Bethe lattice. Note 
that these relations remain valid for the one-dimensional 
chain {K =1). For the triangular Husimi cactus a rela- 
tion similar to JHl is derived in the Appendix. 



B. Green function 

As a first application of the operator identity © we 
present a shortcut to compute the Green function for the 
NN hopping Hamiltonian tiHi on the Bethe lattice, 



G,,(z) = (z|(z-tii/i)-i|j), 



(13) 



with Imz ^ 0. This will fit into Eq. © if we set z 
ti(l + Kx'^)/x, i.e.. 



2ti 



+ - 4:Ktl 



(14) 



with sgn(Im ^ — \Kt\) — sgn(Imz). Since for the 
Bethe lattice there is only one non-self-intersccting path 
connecting sites i and j we have (°L\Hd\2) — Sd,di , whence 



1 



ti 1- 
2K 



E^^ 

d=0 



1 X'^'^+I 

Ul-x'^ 



(K - l)z + Zy/z^ - 4Ktl 



2ti 



■+ 



dij 



(15) 



in agreement with Refs. |2l| and|22. This yields the well- 
known expression for the DOS of Hi, 



pi(A) = --lmGu{X + iO) = , (16) 

where we set \ti\ = l/VT? and p = Z/K = 1 + 1/K. In 
the limit Z ^ oo this leads to the familiar scmielliptic 
density of states, 



lim pi(A) 



2tt 



(17) 



as mentioned in the introduction, with half-bandwidth D 
~ 2. For brevity, here and below we use the convention 
that any square root that appears in a DOS yields zero 
if its argument is negative. 



IV. DENSITY OF STATES FOR 
ARBITRARY-RANGE HOPPING 

The DOS encodes information about hopping param- 
eters and underlying lattice structure. For general tight- 
binding Hamiltonians on the Bethc lattice we now derive 
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an expression for the DOS by using the operator iden- 
tities of the preceding section. In particular for ti-t2 
hopping we expHcitly evaluate the DOS and also the lo- 
cal Green function. The case of exponentially decreasing 
hopping is also discussed. Furthermore we show how to 
construct a tight-binding Hamiltonian from a given DOS. 



Dispersion relations for tiglit-binding 
Hamiltonians 



Since the hopping Hamiltonian Hd on the Bethe lattice 
is a polynomial in Hi [sec Eq. (jl2a|l ]. its eigenfunctions 
are the same as those of Hi, and its eigenvalues can be 
expressed as a function of the eigenvalues of Hi. Ex- 
plicitly, for a general tight-binding Hamiltonian of the 
form (P), i.e., 



(18) 



d=0 



d=0 



we find from Eq. ((na|) that H = T{Hi) with 



ci=0 



K 



Ud{x/2) , 



(19) 



where we used again the scaling td = t'^/K'^/^ and Hd = 
Hd/K'^/^ to facilitate the discussion of the limit Z ^ oo. 
The eigenvalues e oi H are thus related to the eigenvalues 
A of Hi by the "dispersion relation" 



(20) 



which provides a surprising analogy to tight-binding dis- 
persions for crystal lattices with translational symme- 
tries, with A playing the role of crystal momentum. As 
a consequence of Eq. (|20|l the DOS for H can be ob- 
tained by a simple change of variables from that of Hi 
[see Eq. (E 



p(e)= / piiX)Sie~eiX))dX 



(21) 



(22) 



where the sum runs over all solutions Ai(e) of e = e(A). 
Explicit eigenfunctions \6) and eigenvalues X{9) of Hi 
were obtained by Mahan in Ref. 1231. 

We see that using Eqs. l(TH|) - (P^ the DOS can be easily 
obtained for any tight-binding Hamiltonian. At most a 
simple numerical inversion of the polynomial e(A) needs 
to be performed. Below we discuss the case of ti-t2 hop- 
ping as well as an exemplary case of long-range hopping, 
for which this inversion can be performed analytically. 



B. ti-t2 hopping 

The DOS for the Hamiltonian with nearest-neighbor 
(ti) and ncxt-nearest-neighbor (t2) hopping can be found 
from 



Hf t' 



- tiHi + t2-f^2 — 

tlHl + t*iHi -pt;i, 



t%H-2 



(23) 



where again p = Z/ K, and ti ~ t\l\jK and ±2 = ^2/^ 
were scaled appropriately, such that the limit Z — > 00 
can be taken. The dispersion relation and the roots are 
given by 



e{X)^t*2X^ ^t\X-t\v, 
Ai,2(e) 



2t% 



(24) 
(25) 



This yields the DOS 

1 e(tf + 4i^(pt5 + e)) A V4-A,(e)2 



2^ ^tf + 4t*(pt*+e) ^p-A.(e)2/^' 

(26) 

= P-tt,tj(e) = PtI."*2(-£)' (27) 

where Q(x) denotes the step function. In the limit Z 
00 this simplifies to 



e(tf + Mljtl + e)) v/4 - Xf{tY 

-t\±^tf^Ul{t*2+t) 
2t% 



(28) 
(29) 



Note that for 7^ the DOS is asymmetric^ in con- 
trast to the case of random ti-t2 hopping on the Bethe 
lattice.a2& 

We now discuss some limiting cases of Eqs. (|26ll and 
(j2Hll- For t*2^Q (pure NN hopping) we recover Eq. I|16() 
since then one of the roots A^ reduces to e/t\ while the 
other diverges. On the other hand, for tj — > the DOS 
for pure NNN hopping (with 1^2! = 1) is obtained 



Po,i(e) 



\/4 - P - e 



27r (1 - e/Z)VpT7 
In the limit Z — > 00 this reduces to2^ 



Po.-i(-e) . 



Po°?i(e) 



1 V3" 



27r VTT^ 



Po°?-i(-f)- 



(30) 



(31) 



In Fig.[21po,t* (e) (pure NNN hopping) is plotted for sev- 
eral values of Z. This function has a square-root sin- 
gularity on the left-hand side of the band. In the case 
Z ~ 2 the DOS has the same form as for a chain with 
NN hopping, because then the Bethe lattice with NNN 
hopping reduces to two disconnected infinite chains. The 
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FIG. 2: Density of states (001 for pure NNN hopping t2 = 
iiKZ — 1) 7^ on the Bethe lattice, for several numbers of 
nearest neighbors Z; t2 = 1 sets the energy scale. Divergences 
are marked on the horizontal axis. 



DOS for ti-t2 hopping on the triangular Husimi cactus is 
obtained in the Appendix. 

In general, the shape of pt*,t»(e) is determined by two 
dimensionless parameters, the coordination number Z 
and the ratio ^2/^1- view of Eq. H27() it is sufficient 
to consider t^,<2 > 0. In the following we use the pa- 



rameter X = tl/{tl + e [0,1], i.e. 



f* — 
'■2 — 



i* = (1 - x){tl + t^). In Figs. 01 and m we present plots 
of Pti.t'i^) for several values of the parameter x in the 
case Z = 4: and in the limiting case Z — > 00. At small x 
{tl ^ tl) the asymmetry of the DOS develops gradually. 
At some critical value ~ 1/5 a square-root singularity 
appears at the left-hand side of the band. For x < x^, the 
DOS is a smooth function vanishing at both band edges, 
whereas for a; > a;* there is a singularity at the lower 
band edge as well as a cusp within the band with diverg- 
ing first derivative dpt*,t^(e)/de. With increasing x the 
cusp moves continuously to the upper band edge where 
it disappears for x = 1. These features can be under- 
stood with the help of the dispersion relation Eq. (|24|l . 
The solutions Ai,2 [Eq. f^ ] of Eq. contribute only 
if they lie in the interval [—2,2] of nonvanishing pi{X). 
While for x < only Ai contributes, for x > x^ both 
solutions contribute if e lies between the left band edge 
and the cusp, whereas only one solution contributes if e 
lies above the cusp. 

Finally we note that the local Green function can be 
obtained from the DOS by a Kramers-Kronig relation, 
i.e., 




FIG. 3: Density of states 1261 for ti-t2 hopping on the Bethe 
lattice with Z = 4 for selected values oi x = t2/t*, where 
tl,t2 > 0, t* = t'l + t2 sets the energy scale. Divergences are 
marked on the horizontal axis. 



For the case Z 00 this leads to the result 
,G,t.o(^i)-GJ?.o(^2) 



— f 



Zl - Z2 



(33) 
(34) 



where zi,2 i-tl±^/t*^^ + 4t*{z + t*))/2 and G^^o{z) = 
(z — \J z'^ — 4i*2)/(2i*^) is the usual local Green function 
for pure NN hopping as obtained from (|15|l . For pure t2 
hopping reduces to 



1 

2t* 



1 + zlt\ 



(35) 



de. 



(32) 



with the square root given by its principal branch. The 
density of states and Green function for pure NNN hop- 
ping [Eqs. (|31|l and (|35|l ] were previously obtained in 
\ using RPE. 



Ref. 
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X = 0.00 
x = 0.05 
x = 0.10 
x = 0.15 




e/t 



FIG. 4: Same as Fig.El but for Z oo [Eq. 



C. Long-range hopping 

For long-range hopping beyond NNN the expressions 
for A(e) and p{e) become rather comphcated. However, 
the special ease of exponentially decreasing hopping am- 



plitudes, td 



„d-l 



t* , with Hamiltonian 



d=l 



(36) 



allows for an analytical solution (here we assume |w| < 1 
to ensure convergence). The familiar case of pure <i hop- 
ping corresponds to w — > 0. From Eq. (|5Jl it is straight- 
forward to obtain the corresponding dispersion relation 
as 



e«.(A) 



t* 
w 



1 - w^/K 
1 — ivX + w'' 



- 1 



(37) 



1. e., e^(A) = t*\ + 0{w). In the interval -2 < A < 

2, em(A) is a monotonous function with inverse function 
A^(e). The DOS p{e) is again calculated using Eq. H22() . 



t* 1 



27r {wt + t*Y p- X,,{tY/Z 



(38) 




FIG. 5: Density of states for exponentially decreasing long- 
range hopping on the Bethe lattice [Eq. (1361 1. for t* = 1 and 
Z oo. 



where again p = Z / K = \ + 1/K. In Fig. Elp«>(e) is 
plotted for several values of w in the limit Z oo. For 
\w\ < 1 we find a smooth DOS with finite bandwidth, 
while for w — *■ 1 the upper band edge (where Aiu(e) = 2) 
moves to infinity and the DOS decreases like e~^/^ for 
large e. 

It is also interesting to consider only "odd" hopping, 
i.e., between different A and B sublattices, or "even" 
hopping between the same sublattices, although there 
is no immediate physical motivation for this restriction. 
The Hamihonians H' = t* w'*"^^2d-i and H+ = 

t* J^'dLi w'^~^H2d describe odd and even hopping with 
exponentially decreasing amplitudes, respectively. They 
lead to the dispersion relations 



;(A) = r 



(1 - w/K)X 



(1 



w \ (1-1- w)"^ 



w 



A2 



(39) 
(40) 



which are antisymmetric and symmetric in A, respec- 
tively, and yield finite bandwidths except for w ~ I. 
Thus for odd hopping the DOS p~{e) is symmetric as 
in Fig. El whereas for even hopping yO^(e) is asymmet- 
ric (Fig. I?!). We note that /0~(e) is finite for w > 
3 — 2^/2 fa —0.172, while square-root singularities at the 
band edges occur for w < w*. For the latter case an- 
other remarkable feature is that the DOS is nearly con- 
stant in the middle of the band, up to a cusp where it 
rises sharply (see Fig. EJ. On the other hand, pi(e) al- 
ways has a square-root singularity at one band edge [e = 
—pt* /{I + w)], similar to the case of pure NNN hopping. 
Note that for w —1 this singularity moves to infinity 
albeit with vanishing weight; most of the weight is near 
the other band edge in this limit. 

We note that by naively setting w = \/K (without ad- 
dressing questions of convergence) in the above results 
one recovers the case of infinite-range hopping, where 
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FIG. 6: Density of states for exponentially decreasing long- 
range hopping between different sublattices ("odd hopping") 
for t* = 1, Z ^ oo, and w < w,. Divergences are marked on 
the horizontal axis. 
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FIG. 7: Density of states for exponentially decreasing long- 
range hopping between same sublattices ( "even hopping" ) for 
= 1, Z — > oo, and u; < 0. Divergences are marked on the 
horizontal axis. For w = -0.9 the finite peak at the upper 
band edge, as well as the singularity at the lower band edge, 
are outside the plotting range. 



each site is connected to every other site ("complete 
graph" ) and the underlying lattice structure becomes ir- 
relevant. The energy e(A) becomes constant in this case 
and we find /9u,(e) 5{e + t*/VK), in agreement with 
the result for infinite-range hopping on the hypercubic 
lattice, for which the corresponding Hubbard model is 
solvable 4^ 



D. Tight-binding Hamiltonian for arbitrary density 
of states 



In mean-field theories the dependence on the lattice 
typically enters only through the DOS for free particles. 



For simplicity often DOSs with a particularly simple form 
are employed, but when calculating two-particle quanti- 
ties or going beyond mean-field theory it is necessary to 
know the corresponding hopping amplitudes. They are 
also useful for the realistic modelling of strongly corre- 
lated materials where the DOS is obtained from ab initio 
calculations^ In this section we show how to determine a 
tight-binding Hamiltonian H that corresponds to a given 
(one-band) DOS p{e). For hypercubic lattices this task 
has already been addressed in Refs. IssHi^lsoHsiI Here we 
perform such a construction for the Bethe lattice, based 
on the relation l^. 

We start again from the general hopping Hamiltonian 
((T5|l . As noted above, the eigenvalues of H are given by 
the dispersion relation e(A) = T{X) [see Eq. (|19|l ]. We 
can thus use the methods of Bliimeu^^i^ who found a 
similar map between the spectra of H and Hi for the 
hypercubic lattice in the limit of high dimensions. In the 
present case we choose 



.F(A) =Ai(ni(A)), 



(41) 



where ni(A) = J_^pi{\')dX' and ?i(//) = p{e) de. 
Here e^i^ is the lower band edge of p(e) and n{n) denotes 
the inverse function of n(/i). Then > for —2 < 

A < 2, and the DOSs are related by 



pi(A)=.F'(A)p(.F(A)), 



(42) 



which indeed implies Eq. (|21|l . It remains to deter- 
mine the hopping parameters in Eq. H18(l for this choice 
of J-{X). Using the orthogonality of the Chebyshev 
polynomials^ we find that ^^+2/^' ~ ~ ''^d,, where 



(43) 



For the remainder of this section we will consider only 
the limit Z — > cx), for which we have = Ud and ni(A) 
= [1 + \pf{X)]/2 + arcsin(A/2)/7r, and uq = /ep(e) 
=: All yields the first moment of the target DOS. Even 
then the calculation of the hopping amplitudes will typ- 
ically involve numerical integrations. As an analytically 
tractable example we consider the model DOS 



Pie) 



AD 



(44) 



where D is the half-bandwidth and the parameter — 1 < 
a < 1 determines the asymmetry, e.g., the first moment 
is given by Mi = — a/(24Dc^). A constant (rectangular) 
DOS is recovered for a = 0, and a square-root singularity 
at one band edge is present for \a\ = 1. Note that this 
model DOS does not vanish at the band edges, in contrast 
to the DOS proposed in Ref. [51I For the latter, however, 
an analytical calculation of the corresponding hopping 
amplitudes appears infeasible. 
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For the DOS (gH) we have ^(n) = ~D + n/p{-D) 
— 6Min?. In general, upon expanding ^{n) as a power 
series in n, the calculation of t*, involves the coefRcients 



pr{X)U4X/2)[n,{\)rdX 



IT 



sin a; sin(d + l)x 



2x — sin 2x 



2tt 



(45) 



dx , (46) 



where we have used the representation^ Un{x) = sin[(n+ 
1) arccosx]/ sin(arccosx) for the Chebyshev polynomials. 
Evaluating Eq. for s = 0, 1, 2 we obtain for the DOS 

631) 



R{d) 



DR(d) for odd d 
—3MiR{d) for even d ' 

for d = 
for d = 2 
for d = 4 



-1 

3 

35 

-7 
367r2 



(47a) 



(47b) 



-128(d+l) 



;r2(d-2)d2(d+2)2(d+4) 



otherwise 



In Fig. Owe plot the hopping amplitudes vs d, which 
decay slower for large distances than in Eq. due to 
the algebraic behavior of R{d) ^ d"^ . Note that the am- 
plitudes for hopping between the same sublattices (even 
d) are proportional to Mi = 0{a) and thus vanish for the 
special case a = (constant DOS). On the other hand, for 
different sublattices (odd d) the hopping amplitudes are 
independent of the asymmetry parameter a; except for t\ 
« 0.58D they are all negative and drop off quickly (e.g.. 



f* 



14^5 « — 0.033Z?). Regarding tj/t* the quantitative 



difference between hopping for constant and semielliptic 
DOS (pure NN hopping) is thus surprisingly small. 



V. DYNAMICAL MEAN-FIELD THEORY 

We now turn to dynamical mean-field theory 
(DMFT)^5,26,27.28,29,30 fo^, ^j^^, Hubbard model, for which 

the Hamiltonian reads 



-^H — ^^hop + -f^int + ^ext ; 
ija 

Hint = U^ni-^riii , 



(48) 
(49) 

(50) 



where cia and c]^ are the usual fermionic annihilation 
and creation operators for site i and spin a, riicr — cl^c^^ 
is number operator, and iJcxt involves external fields and 
is given below [Eq. H57|l ]. DMFT becomes exact for this 
model in the limit of infinite coordination number, Z 
oom^ Below we will evaluate the DMFT equations for the 
Bethe lattice with arbitrary (in particular ti-t2) hopping. 




FIG. 8: Hopping amplitudes on the Bethe lattice (with Z 
— > oo) corresponding to the DOS 14411 (see inset) for several 
values of the asymmetry parameter a. The "odd" hopping 
amplitudes t*, £3, etc. are independent of a; tg gives the 
center of the mass of the band. 



A. DMFT for the Hubbard model 

We begin with a very brief summary of DMFT for the 
Hubbard model (see Ref. for a review) . Let G and S 
denote the imaginary-time-ordercd Green function and 
self-energy, which arc matrices in site labels i, spin indices 
(T, and imaginary-time slices r (or Matsubara frequencies 
iujn)- They satisfy the Dyson equation 



(51) 



where G*^"-* is the Green function for U = 0. Due to the 
appropriate scaling of tij the self-energy becomes local in 
the limit Z — > 00, 



(S) 



(52) 



and its skeleton expansion depends only on the local 
Green function 

Gian ~ {G)ii^cr,n = ([(G'-"'') — S] ^)ii^a,n ■ (53) 

Therefore Gtan and Sio-n can also be calculated from an 
auxiliary impurity problem^^ with the action 



U / ct^{T)c,^{T)cl(T)c,^(T)dT, (54) 

Jo 



according to 



G^a,l = (c*^(itJ„)c,^(icj„))A 



(55) 
(56) 



so that the so-called Weiss field Gia-n can be eliminated 
and two equations, H53|l and (|55|l . remain for Gian and 
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The DMFT equations |(S2I)-(|SS1) arc thus a closed set of 
equations for the self-energy of the Hubbard model (|48|l . 
and become exact in the limit Z ^ oo. In practice these 
equations are solved by iteration, and the dynamic impu- 
rity problem H54() - (|56|l is solved approximately by numer- 
ical or diagrammatic methods. In addition one must eval- 
uate the self-consistency equation (|53|l . i.e., Gi^n must be 
expressed in terms of Sicr„ and the noninteracting spec- 
trum, preferably involving only the DOS. It is also useful 
to obtain Gian in terms of Gicr„ for use in Eq. ifK^ . Be- 
low these evaluations are performed for ii-t2 hopping on 
the Bethe lattice with Z ^ oo for homogeneous phases 
and phases with broken sublattice symmetry. For ho- 
mogeneous phases S^o-n and Gian are independent of i, 
whereas for broken sublattice symmetry on a bipartite 



lattice one has E,-, 



and Gi, 



Gjan where 7 



= (—1)' = ±1 = A,B depending on the sublattice. We 
add homogeneous and staggered magnetic fields to the 
Hamiltonian, 



(57) 



i.e., with a local field hia = (/if-|-(— l)*/iaf)CT, which allows 
one to detect ferromagnetic or antiferromagnetic phases. 



B. Self-consistency equations 



We now evaluate the self-consistency equation (|53|) for 
arbitrary (scaled) hopping = t*^,jK'^'i/'^ on the Bethe 
lattice. The effective dispersion is then e(A) = with 
DOS /9(e), see Eqs. H19ll - H22() . In particular we consider 
ti-t2 hopping in the limit Z ^ 00 with DOS 



1. Homogeneous phases 

In the homogeneous case all sites are equivalent, as de- 
scribed above. We thus set /laf = but keep ft,f, which 
possibly leads to a ferromagnetic response. The self- 
consistency equation becomes 



Gar, 



■de 



Pi (A) 



(58) 
(59) 



2. Phases with broken sublattice symmetry 

In the case of broken sublattice symmetry the local 
Green function and self-energy depend on i through the 
sublattice index 7. Spontaneous breaking of this sym- 
metry can be detected through the staggered magnetic 
field /laf; We also keep the homogeneous field hi which 
permits phases with Sacto 7^ '^Bsn- 

We start from the eigenbasis of the NN hopping Hamil- 
tonian Hi, i.e., Hi\9) = \{8)\0); these eigenstates may 
later be identified with those of Ref. We use the sub- 
lattice transformation \0) = which yields 
Hi\0) ~ ~~\{0)\0), and introduce wave functions that 
have nonzero amplitudes only on sublattice 7, i.e., \^6) = 
(16*) -|-7|^))/V2. Using the corresponding fermion opera- 
tors c^0^ we now transform the tight-binding part i?hop 
of the Hubbard model with effective dispersion JF(A(6')), 



A(e)>o 




(60) 



where h^^ = {hf+-iKi)a and T±{X) = [J^(A)±J^(-A)]/2. 
Note that the reduced interval for \{9) is analogous to 
the halving of the magnetic Brillouin zone for regular 
crystal lattices. 

For the evaluation of the self-consistency equation H53|) 
we include the self-energy and perform the matrix inver- 
sion. This yields the interacting local Green function as 
a function of the local self-energy 



E 



\{i\ie)\' [z,,n-J'+{m)] 



A(e)>0 



pi{\) [z 

■-fan 



d\, (61) 



where i is any site belonging to sublattice 7 and z. 
iuj„ + ^J■ + h^a — S^CTn- Equation H61|) thus yields the 
DMFT self-consistency equations for arbitrary hopping. 

The general expressions Eqs. (|59|l and (|61|l are the cen- 
tral results of our paper regarding DMFT. Note that their 
derivation involved no counting of lattice paths or other 
combinatorial efforts. Rather they were made possible 
by the dispersion relation ^(A) [Eqs. ((T^ - lf^ ] which is 
due to the operator identity @. 



where z„n ~ iw„-|-/i-(-iT/if — I]„cr, A* is the chemical poten- 
tial, and the dispersion relation ^(A) [see Eqs. ((T^ - (|^ ] 
has been substituted. The general result (|59() thus yields 
the DMFT self-consistency equation in the homogeneous 
case for arbitrary hopping. The special case of ii-t2 hop- 
ping is discussed in Sec. IV GI 



C. t\-t2 hopping 

We now specialize to ti-t2 hopping [with dispersion 
relation (|24|l ] and the limit Z ^ 00. For the homogeneous 
case we find from (|59|) that Gem = G'^ ^, {zan), where the 
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latter function was obtained in Sec. lIVBl [Eas. (|32I35() ]. 

The resulting Weiss field is discussed below [Eq. (|^ ] . 

For broken sublattice symmetry we use Eq. I|61|l with 
T+{X) = (A^ - l)i* and T-{X) = At|. Performing the 
integral in Eq. (|61|) we find 



" 2t* ' 



E 



2(A2 — Af )Ait2^ 



, (62) 



where ±Ai are the poles of the integrand, given by 



Returning to Eq. (|65(l we note that for pure NN hopping 
(^2 = 0) it reduces to the standard cxprcssion^S 



while for pure NNN hopping {t\ — 0) 



2 fcr^ 



(67) 



(68) 



In the former case there is strong sublattice mixing since 
hopping takes place only between different sublattices. 
In the latter case we find no mixing at all since hopping 
is allowed only between sites in the same sublattice. 



A, = ^A±^/A^-B, 

t*2i 



2tf 



S= ( ^ + lU^ + l 



t*2 



(63a) 
(63b) 

(63c) 



and all square roots are given by their principal branches. 
A more compact expression results if one solves for z^fjn 
in terms of Green functions, i.e.. 



(64) 



In the present notation the Weiss field is given by G^^n 
= iijJn + + h^a + G:^an ~ ^jcm- From Eq. follows the 
exact relation between Weiss field and interacting local 
Green function as 



fan 



(1 ^2 X/7' Cr.y'crn)'^ (1 ^2^1(7-0) 



(65) 



Thus for nonrandom ti-t2 hopping the self-consistency 
equation is more complicated than for the two-sublattice 
fully frustrated model with random hopping)S2i^ in 
which all terms in parentheses are absent. Our algebraic 
derivation validates Ref. m where Eq. (^3 was obtained 
using RPE methods, which required the classification of 
many complicated hopping processes. We also note that 
for [/ = Eq. (|?t5|| reduces to the same quartic equation 
that follows from Ref. for the nonintcracting Green 
function [Eq. 

The corresponding equation for the homogeneous case 
is obtained by setting h^i to zero and dropping sublattice 
indices, i.e., 



i*2 



(1 



+ 



)2 (l-f*G,„) 



Gan ■ (66) 



VI. DISCUSSION 

On the Bethe lattice and the triangular Husimi cactus 
the number of possible paths joining two points is de- 
termined only by the topological distance between these 
points. Using this property we derived the operator re- 
lation lO and obtained the spectrum of arbitrary tight- 
binding Hamiltonians, without the need for complicated 
geometrical constructions of other methods. 

In Sec. IIVI the density of states was calculated ana- 
lytically for several classes of tight-binding Hamiltonians 
on the Bethe lattice. As the NNN hopping amplitude 
increases it becomes asymmetric and develops a square- 
root singularity at a band edge. For the Hubbard model 
such a shape of the DOS is known^iiSiS to support 
metallic ferromagnetism away from half-filling, and to 
suppress antiferromagnetism near half-filling. 

In Sec. IIVDI we showed that for any given DOS one 
can determine the corresponding hopping parameters on 
a Bethe lattice. This result is useful in particular in the 
context of dynamical mean-field theory. Namely, with 
this inverse construction it is now possible to mimic any 
kind of van-Hove singularity or other band structure even 
on the Bethe lattice. 

Furthermore, we derived the exact self-consistency 
equations of dynamical mean-field theory in Sec. Ivl for 
arbitrary hopping. In the particular case of ti-t2 hop- 
ping these equations differ from those used in previous 
investigations which are therefore found to ap- 

ply only to random hopping. The exact DMFT setup 
derived in this paper [DOS H28|l and Weiss field H65|) ] can 
be expected to lead to new quantitative results for the 
phase diagram of the ti-t2 Hubbard model on the Bethe 
lattice. This is clear from the exact relation (|65|l between 
the Weiss field and the interacting local Green function, 
since the Green function now appears also in the denom- 
inator, thus leading to resonances in the Weiss field. It 
will be most interesting to study the possible solutions 
obtained within this framework in detail. 

In conclusion, the method and results presented in 
this paper lay the foundation for systematic studies of 
correlated electronic systems on a Bethe lattice. Pos- 
sible future applications include geometrical frustration 
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FIG. 9: Part of the triangular Husimi cactus with Z trian- 
gles connected to each site; here Z = 2. Any two sites are 
connected by a unique shortest path of bonds. The lattice is 
infinite and all sites are equivalent. The shading of sites has 
the same meaning as in Fig. 



in itinerant systems or quantum magnets, and weak- 
localization effects in disordered systems. Indeed, the 
relation lO may also become quite useful for future in- 
vestigations in statistical mechanics. 
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APPENDIX A: THE TRIANGULAR HUSIMI 
CACTUS 

The triangular Husimi cactus is a set of triangles con- 
nected in such a way that each vertex belongs to Z tri- 
angles and each edge belongs to only one triangle (see 
Fig. O. Below we show how to express powers of the 
NN hopping Hamiltonian li\ in terms of the hopping 
Hamiltonian i/d between sites with topological distance 
d [see Eq. (|lb|l ]. For brevity we again set if = Z — 1 and 
p = Z/K, and employ scaled operators Hd = Hd/K'^^'^. 

One can check directly that the Husimi cactus is dis- 
tance regular, i.e., the number of paths between two ver- 
tices i and j depends only on their topological distance. 
This shows that the absence of closed loops as on the 
Bethe lattice is not a necessary condition for this prop- 
erty. The number a)f^ of paths of length n between two 
points with distance dij — d satisfies recursion relations 



very similar to those on the Bethe lattice: for n, d > we 



find a(f+/^ = 2Ka)r'> + a'i'^'' + a^^ a^+i = 2^«" , 
a^rf^ Sod, and ai"^ = 1. Using Eq. Q we thus find for 



(d+2) 



(rf+i) 



-'0 



the first few powers of Hi 

{Hi)^ ^ 2Z1 + Hi + H2 , 
{Hif = 2Z1 + (4Z - \)Hi + 2H2 + i?3 : 
[Hif = (4Z - 1)2Z1 + 5(2Z - l)Hi 
+ (6Z- l)i?2 + 3iJ3 +i?4. 



(Ala) 
(Alb) 

(Ale) 



Note that due to the presence of closed loops there are 
also paths of odd length joining vertices with even dis- 
tance and vice versa. 

We now proceed similar to Sec. IIIII The recursion re- 
lations for the generating function Fd{u) = X]^o '^n^^"" 
can again be solved by an ansata^ Fd{u) ~ f{u) [ug{u)]'^. 
After some algebra we obtain the operator identity 



{l + 2x){l-x) 



1 - x(Ht 



1 



2Kx^ 



(A2) 



d=0 



relating nearest-neighbor and long-range hopping Hamil- 
tonians on the Husimi cactus. (We omit explicit expres- 



sions for the coefficients a"n^ or A)^' .) Note that Eq. (|A2p 
is somewhat more complicated than the corresponding 
relation for the Bethe lattice [Eq. ^] . 

The Green function p3(l can now be obtained directly 
from Eq. lUSJ, similar to Sec lHIBl We find 



(d) 



G,,{z)^{i\{z-tiHi)-'\j) 



2xZti 



(A3) 



where x = [z/ti - 1 + ^/{z/ti - 1)2 - m]/{^K) and 
again denotes the topological distance between sites i and 
j. For Hi (i.e., ti = l/\fK) the DOS is given by 



Pi (A) = 5k 



5{\ + 2) p VS- (A - 1/VX)2 



2^ {p + \/^){2p-\/y/K) ' 

(A4) 



implying a continuous spectrum in the interval 1 / \fK — 
\/8 < A < 1 / \/K + \/8 with an additional delta-peak 
in the case K = 1. Previous derivations of ljA4p can be 
found for K = 1 in Ref. '22 and for A' > 1 in Ref. 
Note that Eq. l|A4p reduces to a semielliptical DOS for 
Z — > oo. 

As an application we now calculate the DOS for the 
Hamiltonian with <i-<2 hopping, i.e.. 



Hti^ti — tHi + t2H2 — i*iHi 

= tHi+t2[Hl - Hi-2Zt] 



(A5) 



due to Eq. ljAla|) . As before we use the scaling ti = 
t\l\fK and t2 — ^jK and find the eigenvalues e of ii in 
terms of the eigenvalues A of iJi, i.e., 



e(A) 



tlx 



+ {t*i~t*2/\/K)X~2pt 



(A6) 
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which has the roots Ai,2(e) = ± / {"^^2) ■ where 

i^^Ul{2pt*+e). (A7) 

Finally the DOS is obtained by changing variables in 
Eq. ||J4)| . 

5{e + 2{tl-tl)) , ^ Pi(A,(e)) 
Pn.tsW +2^^==. (A8) 



The Husimi cactus is not a bipartite lattice, and the sym- 
metry pt*.t;(e) = P-f .t% (e ) holds only for K = 00; this 
can be seen from Eq. We omit plots of l)A8|l since 

its behavior with varying is rather similar to the DOS 
(HHJ for the Bethe lattice. 
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